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> 

^ ' Abstract: We discuss a practical formulation of backward evolution for the CCFM 

small-x evolution equation and show results from its implementation in the new 
^ . Monte Carlo event-generator Cascade. 

p^' 1 Introduction 



In recent years a wealth of experimental results has become available from HERA concerning 
structure functions and final-state properties in deep inelastic collisions (DIS) at small Bjorken 
X and moderate and this has led to interest in theoretical descriptions and predictions of 
the phenomena that are observed. 

DIS at moderate values of x is well described by resummations of leading logarithms of 
transverse momenta (ctglnQ^)", generally referred to as DGLAP physics [ |I], 0, 0, ^. At small 
X we expect leading-logs of longitudinal momenta, (aglnx)", to become equally if not more 
important. However, while the understanding of DGLAP resummations has been mature for 
some time now, despite considerable effort small- a; resummations still remain the subject of 
many theoretical uncertainties and technical difficulties. 

Since many of the measurements at HERA involve complex cuts and multi-particle final 
states, the ideal form for any theoretical description of the data is a Monte Carlo event-generator 
which embodies small-x resummations, in analogy with event generators such as LEPTO [ |5|, 
PYTHIA [ |, HERWIG [ g § and RAPGAP [ | which embody DGLAP resummations. 

In order to build such an event generator two ingredients are required. Firstly one needs 
to know the underlying parton branching equation which, when iterated over many branch- 
ings, reproduces the correct leading logarithms. Secondly one has to find an efficient way of 
implementing the branching equation into a Monte Carlo event generator. 

*Research supported in part by E.U. QCDNET contract FMRX-CT98-0194. 
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As it happens there are several branching equations which are advocated by various groups 
as suitable for describing both inclusive and exclusive properties of small-x DIS. The main 
purpose of this article is to discuss how to formulate one of these equations, the CCFM equation 



[ |T0| |11], |I2|, in a manner suitable for carrying out a backward evolution. This is an 
almost essential requirement if one wishes to efficiently generate unweighted Monte Carlo events 
and modern DGLAP based Monte Carlo generators are always based on backward evolution 
approaches [ |T^ . One of the main results of this paper is that despite the fact that the 

CCFM equation is considerably more complicated than the DGLAP equation, it is possible to 
cast the backward evolution in a form which looks quite similar to the normal DGLAP approach. 
We then show predictions obtained from a new Monte Carlo event generator. Cascade, which 
implements the backward evolution, and compare them with HERA data. 

The paper is structured as follows: in section ^ we discuss briefly the reasons for choosing 
the CCFM equation for the underlying branching. Then in section |^ we discuss the CCFM 
equation itself, and review some details of its implementation in the forward-evolution Monte 
Carlo event-generator Smallx [ ^ |T8[, which has been used to generate the unintegrated 



gluon distribution required by the backward evolution. In section || we discuss the backward 
evolution itself, and then in section |^ show some results. 



2 Why CCFM? 

There are three equations which are commonly used for predictions in small-x DIS: the BFKL 



equation [|T9|, ^, 0, the CCFM equation [|T0|, [11], |T2|, g and the LDC equation [g2[ 0, g, |25 

The BFKL and CCFM approaches are known to reproduce the correct small-x leading 
logarithms for the total cross section. For final-state properties, the derivation of the BFKL 
equation is such that it is not able to guarantee the correctness of the small- a; logarithms. On 
the other hand the derivation of the CCFM equation, based on the principle of colour coherence, 
is such that it guarantees the correct leading logarithms for all final-state observables. (It so 
happens that the CCFM equation also gives a correct description of the final state in the limit 
X 1, another region particularly sensitive to coherence effects, however this is not relevant 
for our purposes). 

This was the situation until a couple of years ago. Recently however it was shown that BFKL 
gives identical leading logarithms of x to CCFM for all final-state observables [ ^ |2S| (the 



physical reasons for this are not fully understood). Since the BFKL equation is quite simple 
(compared to CCFM) one might think that this opens the way to a BFKL-based Monte Carlo 
for DIS. But from the point of view of a sensible description of exclusive quantities, it is not 
just the leading logarithms which matter. For example with the current state of the art, one 
has two options when implementing the BFKL equation in a event generator for DIS. One 
possibility is to use x as the evolution variable — but it turns out that for some observables 
this introduces a weak (but pathological) dependence on one's infrared cutoff in the subleading 
logarithms of Inx)"" (p > 0) [ |2^. The second way of implementing BFKL, which 



resolves this problem, is to use rapidity as the evolution variable. But this turns out to be 
at the expense of introducing subleading logarithms (as In^ Q^)" which violate renormalisation 



2 



group considerations.^ Specifically DGLAP, CCFM and BFKL with x as the evolution variable 
all predict that F2 at small x and large should behave as 



F2ix, Q^) ~ exp InQlnl/ 



'X 



(where = uJJaIt^ and we use fixed for the purposes of illustration). For BFKL with 
rapidity as the evolution variable one has to substitute In 1/a; with the rapidity difference 
between the two ends of the chain Ai] ~ InQ/a;, and this leads to F2 behaving as 



F2{x,Q') ~ exp (^2yaslnQlnl/x + asln^Q 

in contradiction to the DGLAP result. 

The CCFM equation does not suffer from these problems and so forms a good basis for an 
event generator. This does not mean that the CCFM equation embodies all of our knowledge 
about small-x branching, in particular it is known to be incomplete in regions of phase space 
where there is collinear or anti-collinear branching. Furthermore it was discovered in [ |3^ that 
seemingly small modifications of the equation can lead to big differences in its predictions. 
Specifically in a version of the equation without the 1/(1 — z) part of the splitting function, 



the replacement of Qikt — q) with Q{kt — (1 — z)q) in the non-Sudakov form factor, eq. ([12D , 
was instrumental in enabling a fit to F2. However around the same time the exact next-to- 
leading logarithmic (NLL) corrections to small-x evolution, terms as(aslna:)"', became available 
|3^: they state that the power uj governing the growth of quantities like the forward-jet cross 



section at small x should have an expansion a; = 4 In 2q;s — 18.4a;s + . . . (forn/ = 4). On the other 
hand the version of the CCFM equation with the replacement Q{kt — q) ^ G(/ct — (1 — z)q) in 
the non-Sudakov form factor leads to u; = 4 In 2q;s — (75±4)a;g + . . . , which is quite incompatible 



with the known expansion [ Q and so this variant of the CCFM equation can be ruled out. 
The version of the CCFM equation that is studied here has u = 41n2as — (9.2 ± 0.5)a'^ + . . . 
[ ^ — i.e. its NLL corrections are somewhat smaller than the true NLL corrections, however 
they are at least of the right order of magnitude. 



There exists also a third evolution equation, the Linked Dipole Chain (LDC) [|22|, |23|, |2^, |25 
whose characteristic is that some of the initial state radiation has been moved into the final 
state. This has two consequences — on one hand the simplifications that ensue mean that it is 
quite easy to correctly implement a symmetry between branching up and down in transverse 
scale (or equivalently one obtains the same predictions whether one evolves from the virtual 
photon or from the proton). Such a symmetry leads to the implicit inclusion of important NLL 
terms associated with large transverse logs in the anti-collinear limit (terms Q;g/(1 — 7)'^ where 
7 is the Mellin variable conjugate to squared transverse momentum). These terms are entirely 
missing from standard LL BFKL, while in CCFM they are present, but with half the correct 
coefficient.0 The second consequence of the manner in which initial state radiation has been 
moved into the final state, is that the LDC has slightly different small-x leading logs compared 
to BFKL {lu ~ 3.23as + • • • for LDC as compared to ~ 2.77a^ + ... for BFKL and CCFM 



^For certain restricted applications in hadron-hadron scattering, where BFKL generators have been advo- 
cated and developed [ ^ |3^, these double- logarithms may not be too important because of the nature of 
the observables studied. 

^There have been attempts to implement the so-called kinematic constraint |3^, which in BFKL does 

lead to the correct NLL behaviour for 7 — > 1. But in CCFM the situation is more subtle and a straightforward 
inclusion of the kinematic constraint still gives the wrong 7 — > 1 limit. 
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Figure 1: Kinematic variables for multi-gluon emission. The t-channel gluon four-vectors are 
given by ki and the gluons emitted in the initial state cascade have four-vectors pi . The upper 
angle for any emission is obtained from the quark box, as indicated with S. 



— it is possible though to modify LDC so that it has leading logs which are much closer to 
the BFKL ones [ 0). Therefore LDC has both advantages and disadvantages compared to 
CCFM. The question of which matters more will almost certainly depend on the nature of the 
observable under study. 

So there remains much room for further theoretical progress in the description of small-x 
evolution. However given the current situation, one of the more viable options is the CCFM 
equation, hence our interest in implementing a practical CCFM-based event generator. 



3 The CCFM evolution equation 



The implementation of CCFM [ |T0], |TT], |T^, |13| parton evolution in the forward evolution Monte 
Carlo program Smallx is described in detail in [ |l^ Here we only concentrate on the 



basic ideas and discuss the new treatment of the non-Sudakov form factor 
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Figure |l| shows the pattern of QCD initial-state radiation in a small-x DIS event, together 
with labels for the kinematics. According to the CCFM evolution equation, the emission of 
partons during the initial cascade is only allowed in an angular-ordered region of phase space. 
The maximum allowed angle S is defined by the hard scattering quark boxj^ which connects 
the gluon to the virtual photon. In terms of Sudakov variables the quark pair momentum is 

■^Strictly there are two maximum angles, corresponding to the directions of the quark and the anti-quark, 
and which each hmit roughly one half of the radiation. 
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written as: 



Pq+Pq = Y{Pp + Spe) + Qt (l) 

where Pp and Pe are the proton and electron momenta, respectively and Qt is the transverse 
momentum of the quark pair. Similarly, the momenta pi of the gluons emitted during the initial 
state cascade are given by (here treated massless): 

Pi = ViiPp + iiPe) +Pti , 6 = (2) 

with Ui = {1 — Zi)xi-i and Xi — ZiXi-i and s = {pp+PeY being the total electron proton center 
of mass energy. The variable is connected to the angle of the emitted gluon with respect to 
the incoming proton and ,Xj and Ui are the momentum fractions of the exchanged and emitted 
gluons, while Zi is the momentum fraction in the branching — and pu is the transverse 

momentum of the emitted gluon. 

The angular-ordered region is then specified by: 

6 < < • • • < < S (3) 

which becomes: 

Zi-iqti-i < qti (4) 
where we use the rescaled transverse momenta qu of the emitted gluons defined by: 

Qti ^ Xi^iy/s^i ^ zr^^ (5) 

i Zi 

In Smallx, the initial state gluon cascade is generated in a forward evolution approach 
from an initial distribution of the kt unintegrated gluon distribution according to: 

xAo{x, kl) = N-bhl-xY- exp {-kl/kl) (6) 

where N is a normalization constant. The exponential distribution in k^ has a width which will 
be set to /cq = 1 GeV^. The input gluon distribution needs to be adjusted to fit existing data, 
but it turns out that the small x behaviour of the structure function F2 is rather insensitive to 
the actual choice of xAo{x,k'^Q) and only the normalization N acts as a free parameter, with 
the constraint: 

J xAo{x, k^Q)dxdk^Q — J xGo{x,Q'^)dxc^0.5 (7) 



which gives: 



J xAo{x, k'^o)dxdk^Q = N-5- J {1 - xfdx = N (8) 



The initial state branching of a gluon ki into another virtual (t-channel) gluon /cj+i and a 
final gluon pj+i (treated on mass shell) is generated iteratively from the initial gluon distribution 
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at a starting scale Qo- The probabihty for successive branchings to occur is given by the CCFM 
sphtting function [ 0, |ll], |T|, 0] : 

dPi = PUzi, q^, kl) ■ Asdzi — ^ ■ 6(gj - Ziqi^i) ■ 0(1 - - e^) (9) 

vrg. 

with qi = pti/{^ — Zi) being the rescaled transverse momentum of the emitted gluon %. The 
fractional energy of the exchanged gluon % is given by Xi and the energy transfer between the 
exchanged gluons i — \ and i is given by Zi = Xi/xi-i. A coUinear cutoff = Qo/qi is introduced 
to regularize the — z) singularity. The Sudakov form factor is given by: 

A,{qi,Ziq,^i) = exp - / / dz^^ ^ 10 

with as = = For inclusive quantities at leading-logarithmic order the Sudakov form 
factor cancels against the 1/(1 — z) collinear singularity of the splitting function. Coherence 
effects are taken into account by angular ordering qi > given by the first G function in 

eq.(|]). The cascade continues until qi reaches the limiting angle defined by g = set by 

the partons from the hard scattering matrix element. 

The gluon splitting function is given by:0 

Qs( g- (1 - Z if) , as{kl) ^ 

I - Zi 

where the non-Sudakov form factor A^^ is defined as: 

log A„, = -a,(4) /' ^ / ^e{ku - q)Q{q - z'qu) (12) 
Jo ^ J Q 

The principle difference compared to the corresponding DGLAP splitting function is the ap- 
pearance of the non-Sudakov form factor A^^, which screens the 1/z singularity in eq.([TT|). It 
can be expressed as [ PB| : 

^0 ^ , _ „^ f kli 



Pi = "-^^r: + ^^nsiz., ql kl) (11) 



log A„, = -a,{ki;) log ^ log ^ (13) 

\ZiJ \ZQZiqi J 

where 

( 1 if ku/qi > 1 

^0="! ku/qi if Zi< ku/qi < I 
[ Zi if ku/qi < Zi 

which means that in the region ku/qi < Zi we have A„s = 1, giving no suppression at all.| 

The CCFM equation for the unintegrated gluon density can be written [ |3^ as an 

integral equation: 

/dz f d^Q ~ fx \ 

— / —\Q{q-zq)As{q,zq)P{z,qX)-A[-X.(l) (14) 
z J T^q ^z ' 



''Actually the 'correct' scale for ttg, as suggested by the NLL corrections to BFKL, is probably — z^)^ in 
both terms, with a corresponding modification of the non-Sudakov form factor. However for simplicity, at this 
stage we have retained the scale for that was present in the original formulation of Smallx. 

^We note that in the original version of Smallx a simplified version of the non-Sudakov form factor was 
used, which however did not give the right answer for kt < q. 
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with k'f. = \kt + {1 — z)q\ and with q being the upper scale for the last angle of the emission: 
q > Znqn, In > ^n-iqu-i, <?! > Qo- As before q is used as a shorthand notation for the 
2-dimensional vector of the rescaled transverse momentum q = qi = pt/{l — z). The splitting 
function P{z,q,kt) is defined in eq.(|lT]) and the Sudakov form factor As{q,zq) is given in 
eq.(0). 

In [ |13| a differential form of the CCFM evolution equation is given, which is obtained from 
eq. ( p!^ by dividing both sides by As{q, Qq) and then differentiating with respect to q: 

2 d xA{x,kt,q) f d(f) P{z,q/z,kt) 



with x' = x/z and k'^. = {1 — z)/zq + kt and where g'has an azimuthal angle 0. In deriving this 
equation we have exploited the fact that the Sudakov form factor can be written as 

^s{q,zq) = — — -. (16) 

As{zq,Qo) 

For (|15|) (and the backward evolution formalism which follows from it) to be correct Ao{x, kt, q), 
must be of the form 

Ao{x, kt, q) = Ao{x, kt) As{q, Qo) ■ (17) 



4 Backward evolution: CCFM and CASCADE 

The forward evolution procedure as implemented in Smallx is a direct way of solving the 
CCFM evolution equation including the correct treatment of the kinematics in each branching. 
However the forward evolution is rather time consuming, since in each branching a weight 
factor is associated, and only after the initial state cascade has been generated completely can 
it be decided whether the kinematics allow the generation of the hard scattering process. Quite 
often, a complete event has to be rejected. 

A more efficient procedure to adopt in a full hadron-level Monte Carlo generator is a back- 
ward evolution scheme, analogous to that used in standard Monte Carlo programs [ [l^, |^ 



using a DGLAP type parton cascade. The idea is to first generate the hard scattering process 
with the initial parton momenta distributed according to the parton distribution functions. 
This involves in general only a fixed number of degrees of freedom, and the hard scattering 
process can be generated quite efficiently. Then, the initial state cascade is generated by going 
backwards from the hard scattering process towards the beam particles. In a DGLAP type 
cascade the evolution (ordering) is done usually in the virtualities of the exchanged t-channel 
partons. 

According to the CCFM equation the probability of finding a gluon in the proton depends 
on three variables, the momentum fraction x, the transverse momentum squared k^ of the 
exchanged gluons and the maximum angle allowed for any emission q = XnVsE. To calculate 
this probability, in addition to the details of the splitting, eq. (P) one needs to know the 
unintegrated gluon distribution A{x, k, q) which has to be determined beforehand. 

Given this distribution, the generation of a full hadronic event has three steps, implemented 
in a new hadron-level Monte Carlo program. Cascade: 
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Firstly, the hard scattering process is generated, 

0- = y dk'fdXgA{Xg, kl, q)a{'j*g* qq) , (18) 



using the off-shell matrix elements given in [ ^ p. 178 ff], with the gluon momentum (in 
Sudakov representation): 

k = XgPp + XgPe + kt^ XgPp + ^ . (l9) 

The gluon virtuality is then — A;^ ~ k'j. 

The initial state cascade is generated according to CCFM in a backward evolution ap- 
proach (described in the next section). 

The hadronisation is performed using the Lund string fragmentation implemented in 
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Strictly speaking, before the last step one should also include angular-ordered final-state 
radiation from the intial-state gluons. For simplicity, in this 'proof of concept' version of the 
generator, this is currently left out. 

In the backward evolution there is one difficulty: The gluon virtuality enters in the hard 
scattering process and also influences the kinematics of the produced quarks and therefore the 
maximum angle allowed for any further emission in the initial state cascade. This virtuality 
is only known after the whole cascade has been generated, since it depends on the history of 
the gluon evolution. In the evolution equations itself it does not enter, since there only the 
longitudinal energy fractions Zi and the transverse momenta are involved. This problem can 
only approximately be overcome by using k'^ = k^/{l — Xg) for the virtuality which is correct 
in the case of no further gluon emission in the initial state. 

The Monte Carlo program Cascade can be used to generate unweighted full hadron-level 
events, including initial-state parton evolution according to the CCFM equation and the off- 
shell matrix elements for the hard scattering process. It is suitable both for photo-production 
of heavy quarks as well as for deep inelastic scattering. The typical time needed to generate 
one event is ~ 0.03 sec, which is similar to the time needed by standard Monte Carlo event 
generators such as LEPTO [ || or PYTHIA [||. 



4.1 The unintegrated gluon density 

The unintegrated gluon density xA{x, k^, q) is obtained from a forward evolution procedure as 
implemented in Smallx [ |TH[. Due to the complicated structure of the CCFM equation, no 
attempt is made to parameterize the unintegrated gluon density. Instead, the gluon density is 
calculated on a grid in logx, logfcj and logg of 50 x 50 x 50 points and then linear interpolation 
is used to obtain the gluon density at values in between the grid points. 

From the initial gluon distribution as used in Smallx (including the same collinear cutoff 
and normalization) a set of values xq i and kto i are obtained by evolving up to a given scale log q 
using the forward evolution procedure of Smallx. This is repeated 10'' times thus obtaining 
a distribution of the unintegrated gluon density XnA{xn,k^^,qtn) for the slice of phase space 
with a given q {q > qtn)- To obtain a distribution in logg, the above procedure is repeated 
from the beginning 50 times for the different grid points in log g up to g = 1800 GeV. 
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4.2 Backward evolution formalism 



In the backward evolution we start from the quark box, with an upper angle given by H and 
a gluon four-vector kn (see Fig. |1]) and go successively down in the ladder until we end up at 
gluon ko. Thus the first step is to reconstruct from S and fc„ the vectors = Pn/(1 ~ Zn) and 
kn-i = ki, with Zn = Xi/xn- lu the next step ki, g„ play the role of kn, H from the first step. In 
the further steps fcj-i, qi play the role of ki, g^+i and so on, until the gluon ko is reached. 

The differential form of the evolution equation eq.([T5|) gives the (non-normalized) proba- 
bility [ , that during a small decrease of q, a t-channel gluon k' with momentum fraction x' 



becomes resolved into a t-channel gluon k with momentum fraction x = zx' and an emitted 
gluon qi. During a small decrease of q, a gluon k may be unresolved into a gluon k'. The 
normalized probability for this to happen is given by 

AMQo) JxA{x,k^q)\ ^df r ^^dl^^^^_^^^^^^x'A{x',K,q/z) ^20) 



xA{x,kt,q) V ^s{q,Qo) J q^ J 27: ' ' xA{x,kt,q) 
This equation can be integrated between q and q giving: 

[ A{x,kt,q) As{q,Qo)\ H dq'^ f a h( '/ , .x' A{x' , k'^, q' / z) 

log u .^ A ^„ ^ ^ = - / -JT dz— P{z, q / z, h) ^ . , ,^ (21) 



A{x,kt,q) l\s{q,Qo) ) g' 7 27r ' ' xA{x,kt,q' 

Thus the probability for no radiation in the angular ordered region between q and q is just 
given by a new effective form factor for the backward evolution: 
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-n (- \ I I dq f , d(j) ~ xA{x,k^,q z)\ 

Vnoradiq, q) = exp [ - dz— P{z, q z, kt) — — (22) 

V Jq 9 J 27r xA[x,kt,q') J 

In principle one could equally well just use 

A{x,kt,q) /^s{q,Qo) 
A[x, kt,q) As{q,Qo) 

which is more akin to the backward evolution approach of [ |T6|. One can see that this is the 
correct probability since from eq. (p!4D it corresponds to the fraction of A{x, kt, q) which comes 
from angles below q. Though eq. ( [2^ ) looks more complicated than eq. (p3D, it turns out that 
the former is numerically more suited to our particular situation, because it is less sensitive to 
imprecisions and irregularities of A{x, k, q) (which we recall is generated by a forward evolution 
Monte Carlo approach). The standard DGLAP backward evolution equation would be obtained 
from eq. ( ^Of ) by setting A„s = 1 in the splitting function P and replacing the argument q/ z in 
the parton density function in the r.h.s. of eq. (^) with q. 

In the CCFM backward evolution, starting with the gluon kn, we need to reconstruct the 
momentum of the next emitted gluon g„ as well as that of the next exchanged gluon, fcn-i- 
This is done as follows. We start with g = S. Then 

a. We check whether there is any evolution to be done. The probability that the evolution 
should stop straight away is given by 

A(a^n,fe?n,g) ^24) 

A{Xn, kin, Q) 

i.e. the fraction of the unintegrated gluon distribution that comes from the initial distri- 
bution. 
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b. If the evolution continues then the quantity q' = \znqn\ is determined by choosing a 
random number R uniformally distributed between and 1 and solving VnoradiQ, <?') = R 
for q'. 

c. The values of Zn and (pn are then chosen randomly according to the distribution given in 
the inner integral of eq. (pSf). In doing so we implicitly have to reconstruct kn-i as well. 



This completes the reconstruction of a single branching. The procedure is then repeated with 
q = q-n and n ^ n — 1, and so on, until the evolution stops (step a). 



In practice the numerical calculation of the integrals in the effective form factor eq.(22) 
would be too time consuming, since they involve the non-Sudakov form factor in the splitting 
function and also the ratio of the structure functions. A simple solution to the problem is the 



veto algorithm described in [ |15| . The essential point is to find a simple analytically integrable 



function, which is always larger than the integrand in eq . (|2^) . The splitting function P{z, g', kt) 
can be replaced with: 

PT{^^, kl) = + ^i^^ > P,, (25) 



The limits on Zi are given by: 



CC 2, * 



which is a larger range than the true one: Zi < 1 — Qo/qi (but qi is not determined at this stage). 
Next, the structure functions which appear in the Sudakov form factor are replaced with their 
maximum and minimum values for x and x' > x. Thus a simple analytically calculable form of 
the Sudakov form factor is obtained: 



After qi and Zi is generated, the true limits on z can be applied: Xi < Zi < 1 — Qo/qi. If a lies 
outside the true region, a new set of q^ and Zi is generated. Having generated the branching 
variables according to Zt, ku-i and qi according to eq.(P^) and eq. (pSj) , a branching is accepted 
with a probability according to the ratio of the integrands of V (via eq.(|^)) and V^^/J^^ (via 
eq.(P^)), as formulated in the veto algorithm [ P^. 

As mentioned already above, the true virtuality of the t-channel gluons can only be recon- 
structed after the full cascade has been generated. By going from the last gluon (closest to the 
proton), which has virtuality /cq = kto/i^ ~ ^o); forward in the cascade to the hard scattering 
process, the true virtualities of the kf are reconstructed. At the end, the gluon entering to 
the quark box will have a larger virtuality than without initial state cascade. Thus a check is 
performed as to whether the production of the quarks is still kinematically allowed. If not, the 
whole cascade is rejected, and the event without the cascade is kept. This typically happens 
about 1% of the time. 
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4.3 Comparison with the forward evolution cascade in Smallx 

In this section, the reconstruction of the parton level cascade obtained in the backward evolution 
approach, as described above, is compared to that of the forward evolution in Smallx. At 
parton level, both approaches are expected to be identical, but small differences can occur due 
to the finite grid size used to define of the unintegrated gluon density. As already mentioned 
the virtuality of the gluon entering the hard scattering process is only known after the complete 
reconstruction of the initial state cascade, which could also result in small differences to the 
forward evolution approach. In the following comparison, we have used bb photoproduction 
at y/s = 300 GeV. In Fig. ^ the cross section as a function of the rapidity (all rapidities 




Figure 2: Comparison of the cross section obtained from the backward evolution Monte Carlo 
Cascade (solid line) with Smallx (dashed line) both at parton level only. The upper plots 
show the cross section as a function of the quark rapidity rjg, the quark transverse momentum 
Pt and the transverse momentum of the quark pair p^"'^^ . The lower plots show the cross section 
as a function of the gluon transverse momentum kt, and the multiplicity and transverse energy 
flow of the gluons from the initial state cascade as a function of the rapidity rj. 

are given in the laboratory frame) of the quarks rjg, the quark transverse momentum pj, the 
transverse momentum of the quark pair p^"'^^ and the gluon transverse momentum kt obtained 
from Cascade (solid line) are compared to the ones obtained from Smallx (dashed line). 
The quantities related to the hard scattering matrix element agree very well. Also shown in 
Fig. is a comparison of the multiplicity and the transverse energy flow as a function of rapidity. 
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and perfect agreement is again found between the backward and forward evolution approaches. 
In Fig. ^ a more detailed comparison of the kinematics in the initial state cascade is performed. 




Figure 3: Comparison of quantities of the initial cascade obtained from the backward evolution 
Monte Carlo Cascade (solid line) with Smallx (dashed line) both at parton level only; (a) 
shows the splitting variable z, (b) gives the transverse momentum kt, (c) shows the transverse 
momentum of the emitted gluon qt, and (d) shows the ratio kt/qt- 

We compare the values of the splitting variable z, the transverse momenta kt and qt as well as 
kt/qt within the two approaches. 

Thus we have described a backward evolution approach, which is fast and already imple- 
mented in the hadron-level Monte Carlo program Cascade, and which reproduces perfectly 
the parton level configurations obtained from Smallx. This is the first time that a practical 
CCFM-based small-x Monte Carlo event generator has been constructed. 



5 Results 

In this section we compare predictions from CASCADE with recent measurements made at 
HERA. The free parameters of the initial gluon distribution were fitted to describe the structure 
function F2(x, Q^) in the range x < 10~^ and > 5 GeV^. The structure function F2(x, Q^) 
as calculated from CASCADE is compared to a larger range of data in Fig. ^. Reasonable 
agreement with the data is observed at small x and Q^. Deviations from the data are seen at 
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larger x and Q^, which can be attributed to the quark contributions, which are still missing 
here. 

In Fig. ^ the cross section predicted for forward-jet production is shown and compared to 
measurements done at HERA. We observe a reasonable description of the data. However, the 
prediction lies above the data, which could indicate that some relevant sub-leading effects are 
still missing. 

Studies based on a variety of QCD-based Monte Carlos have demonstrated that the high 
Pt tail of charged particle transverse momentum spectra is sensitive to small x dynamics of 
the parton radiation and that there is relatively little dependence on the particular model of 
hadronisation that is used [ Fig. |^ shows the distributions of charged particles as 



measured by HI [ ^ for DIS events with 3 GeV^ < < 70 GeV^ in the rapidity range 
0.5 < 7] < 1.5. The prediction from CASCADE gives a good description of the data, whereas 
DGLAP based Monte Carlo calculations fail to describe the data at small x and large Pf 

Effects of small x parton dynamics could also show up in photo-production of charm mesons. 
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Figure 5: The cross section for forward-jet production obtained from the Monte Carlo Cas- 
cade at hadron level (solid line); (a — c) The cross section for forward-jet production as a 
function of x, for different cuts in pt compared to HI data [\4^J (a — b) and compared to ZEUS 
data (c); (d) The cross section for forward-jet production as a function of E^/Q"^ compared 

to [dw. 



Recent measurements [ show significant deviations from LO and NLO QCD calculations 



and also from hadron-level Monte Carlo predictions. In Fig. ^ we show the cross section of D* 
production as a function of the transverse momentum p^* using CASCADE and compare it with 
the measurement of ZEUS [49]. For comparison we also show the prediction from Rapgap. 



In Fig. 1^6 the cross section is shown. In Fig. |^ we show the D* cross section as a function of 
the pseudo-rapidity rj^* for different regions in pt. In all distributions of D* photo-production 
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Figure 6: The transverse momentum distribution of particles in different bins of rapidity. The 
prediction from CASCADE (solid line) at hadron level is compared to the measurement of HI [ 
For comparison the prediction from the DGLAP based Monte Carlo Rapgap (dashed line) 
is also shown. 



we observe a good description of the experimental data, whereas there is disagreement between 
the data and the Rapgap Monte Carlo predictions. 
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Figure 7: The differential cross sections of D* photo-production [ for da/dp^* (a) and 
da/dx^ (b). The solid line shows the prediction from Cascade while the dashed line shows the 
prediction from Rapgap. 



5.2+_li nb. 



Also interesting is the cross section for bb production as measured at HERA by HI [ |5D| : 
a{ep e'bbX) = 7.1 ± 0.6{stat.)tli{syst.) nb 
to be compared with the prediction from Cascade: 

a{ep e'bbX) 

The central value is given for = 4.75 GeV, and the errors are those associated with a 
variation of mj, by =p0.25 GeV. The NLO calculations predict a cross section which is about a 
factor of 2 below the measurements. The ZEUS collaboration published a measurement of bb 
production [jSTjl for the region pt > 5 GeV and \rii,\ < 2: 

a{ep e'bbX) = 1.6 ± OA{stat.)toi{syst.)toi{ext.) = l.Qto.rt nb 

where we have added all errors quadratically in the last expression. This can be compared with 
the prediction from CASCADE: 

a{ep e'bbX) = 0.88 ± 0.08 nb. 

Again we have used m?, = 4.75 =F 0.25 GeV. The predicted cross section is still within the total 
error of the measurement. The NLO prediction as given in [ |^ is cr = 0.641°;}^ nb, where 
the errors are those associated with a variation of = 4.75 GeV by =f0.25 GeV and of the 
factorisation and renormalisation scales by a factor of 1/2 and 2. 



6 Conclusion 



We have shown that a backward evolution approach using the CCFM evolution equation is 
possible, and that it works, producing the same results as the forward evolution used to solve 
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Figure 8: The differential cross section of D* photo-production [ \4^J da/drj^ for different 
regions of p^ . The solid line shows the prediction from Cascade and the dashed line is the 
prediction from Rapgap. 

the CCFM equation. The advantage of backward evolution is the easy implementation into a 
full hadron-level Monte Carlo program Cascade, which is compared to recent measurements 
of small X hadronic final state properties at HERA. We have found that all small x signatures 
can be reasonably well described within one consistent approach. By performing only a fit to 
the structure function F2{x,Q'^) we obtain simultaneously a good description of a variety of 
processes which could not be described within DGLAP: the forward-jet cross section, the high 
Pt particle spectra, charm production and a reasonable agreement with the measured bb cross 
section. This shows that the CCFM evolution equation is indeed a good starting point for a 
consistent description of small x phenomena. 
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